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Abstract 

We perform a global analysis in the framework of two active neutrino os- 
cillations of all solar neutrino data, including the recent SNO day and night 
spectra (comprised of the charged current (CC), elastic scattering (ES) and 
neutral current (NC) events), the Super-Kamiokande (SK) day and night 
spectra (from 1496 days) and the updated SAGE results. We find that the 
Large Mixing Angle (LMA) solution is selected at the 99% C.L.; the best-fit 
parameters are Am? = 5.6 x 10~^ eV^ and 6 = 32°. No solutions with 9 > 7r/4 
are allowed at the 5a C.L. Oscillations to a pure sterile state are excluded at 
5.3(7, but a sizeable sterile neutrino component could still be present in the 
solar flux. 
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The end of the era of the solar neutrino problem is on the horizon. Neutrino oscillations 
are a compelling explanation of the solar neutrino deficit relative to the Standard Solar 
Model found by several experiments |0-§1. The SNO experiment has provided cogent 
evidence in favor of this hypothesis by separately measuring the incident z/g flux via the CC 
reaction and the total active neutrino flux via the NC reaction in the same energy range [0 . 
With the continued accumulation of solar neutrino data, the LMA solution has become 
the most favored. If LMA is in fact the solution, the ongoing KamLAND experiment 



will provide a precise and SSM-independent measurement of the oscillation parameters by 
measuring the suppression and distortions of the z/g flux emerging from the surrounding 
nuclear reactors |[Tl|]. Aside from a possible complication of partial oscillations to sterile 
neutrinos, the thirty year old solar neutrino problem will be solved. 

Here we assess the extent to which the recent SNO results home in on an unique oscil- 
lation solution by performing a global analysis of all solar neutrino data in a two neutrino 
framework. There is little difference between a three neutrino analysis and an effective two 
neutrino analysis because the mixing between the solar and atmospheric scales is known 
to be small |]12[. We first show through a model- independent analysis that pure active to 



sterile oscillations are excluded at high confidence and then proceed with a global analysis 
involving only active flavors. 

A universal sum rule that holds for arbitrary active-sterile oscillation admixtures has 



been derived [|T3l which relates the NC flux at SNO to the CC flux at SNO and the neutrino 



flux measured at SK with elastic scattering: 

<l>NC = ['^SK-{l-r)<!>cc]/r, (1) 

where r = ajj^^^v^/ (Ty^ is the ratio of the u^^r to z/g elastic scattering cross sections on electrons. 
Above the 5 MeV threshold of the two experiments, r = 1/6.48 and the relation becomes 

= 6.48<l>5i^ - 5.48<l>cc ■ (2) 



The values from the data (assuming an undistorted B spectrum 10-0), for the left and 
right-hand sides of Eq. (0) are, respectively, 

5.09 ±0.62 , 5.58 ±0.71 (3) 

showing agreement between the data and the sum rule. The fact that the two values have 
similar percentage uncertainties (~ 12-13%) shows that the accuracy of the SNO NC mea- 
surement is already comparable to that which can be inferred from the SNO CC and SK 
data. 

The SNO NC rate is a measure of the flux of active neutrinos in the high energy part 
of the solar neutrino spectrum. If an active-sterile admixture is responsible for the solar 
neutrino deficit. 



^SSM ^SSM 

with measured value (see Table D, 

^NC 



1.01 ±0.12. (5) 

^SSM 
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TABLE I. Solar neutrino data relevant to our analysis. Note that we use the SNO day and 
night spectra in place of the SNO CC, NC and Aijjy values and SK day and night spectra instead 
of the SK rate. The SNO CC, NC and Adn numbers were extracted from the SNO day and night 
spectra assuming an undistorted spectrum and are given here only for reference as is the value 
of ^ssM- 1 SNU is 1 interaction/s/10^^ atoms of the neutrino absorbing isotope. 



Measurement 


Homestake 


2.56 ±0.23 SNU 






GALLEX+GNO 


73.3 ± 4.7 ± 4.0 SNU 






SAGE 


70.8 ± 5.3 ± 3.5 SNU 






SK 


2.35 ± 0.02 ± 0.06 x 10^ 


cm" 


-2 S-^ 


SNO CC 


1.76 ±0.06 ±0.09 X 10^ 


cm" 


-2 


SNO NC 


5.09 ± 0.44 ± 0.45 x 10^ 


cm" 


-2 S-1 


SNO Adn 


0.07 ±0.049 ±0.013 






^SSM 


5.05(1 ±0.18) x 10*^ cm 


-2 s- 


-1 



Here /? is a normalization of the flux with respect to the central value of the SSM 
prediction ^ssm = 5.05 x 10^ cm~2s~^ and (3^us is the total sterile neutrino flux to which the 
electron neutrinos oscillate^. Since the measured ^nc is consistent with the SSM prediction, 
we impose the SSM constraint /3 = 1 ± 0.18, and obtain 

P{ve Vs) = -0.01 ± 0.22 . (6) 



SSM 



While the central value of the SNO NC rate suggests a solution with oscillations only to 
active flavors, the uncertainty is too large to rule out a substantial sterile neutrino flux on 
this basis aloneQ. 



To see this freedom from another angle, we return to the approach of Ref. [|1^]. Of the 
neutrinos that oscillate, the fraction that oscillate to active neutrinos is 

sm a = — — — — = 1.01 ± 0.34 . (7) 

P^SSM — '^CC 

Thus, the evidence for transitions to active flavors is at the 3a CL. However, large sterile 
fractions are allowed even at the 2a CL. See Fig. ^ for illustration. The dark-shaded 
and light-shaded regions enclose the values of ^nc and allowed by the SSM at la 
for sin2 a = 1 and sin2a = 0.5, respectively. The lines through the center of these bands 



^As in the SNO analysis we conservatively adopt the SSM ^B flux of BPB2000 since the 
S'i7(0) determination used in Ref. [|l^] is being reanalyzed [ll6| . 

2Assuming an LMA solution and a small sterile admixture, a combination of KamLAND and 
SNO CC data can determine /3 to 11% [|l^. This uncertainty is still too large to eliminate the 
possibility of a significant sterile component in the solar flux. 
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correspond to the central value of the SSM flux prediction. The region above the diagonal, 
$7vc7 = $CC) is forbidden because ^cc > ^nc impossible. The measured SNO CC and 
NC rates are marked by a cross. Even a doubling of the widths of the SSM bands and the 
SNO rate uncertainties (effectively 2a) allows large sterile fractions. 

We emphasize that a large sterile neutrino flux is viable not only for an analysis of 
the total rates, but also when the SK day and night spectra are included (contrary to the 
assertion in Ref. ||TB|]). For the issue in question, the effect of imposing the SSM flux 
constraint is equivalent to including the day and night spectra: large values of f3 are not 
allowed and the best-flt value of sin^ a is close to pure active oscillations. That this is the 



case can be seen from Fig. 4 of Ref. [|T3| ( B flux constraint applied) and Fig. 2 of Ref. |T7 
(day and night spectra used). The reason for this correlation is that the day and night 
spectra rule out a large day-night effect, which is the same region of the LMA solution (low 
Am^) that favors large P for small sin^ a. Therefore in our rate analysis, imposing the ^B 
flux constraint should yield the same effect as including the SK day and night spectra, and 
we see from Eq. (^ that a large sterile flux is still allowed. 

In the above analysis, pure active and pure sterile oscillation solutions are treated on an 
equal footing. If instead we use the a priori /3-independent criterion that for a pure sterile 
oscillation solution, 

"^NC = $CC7 (8) 

(the diagonal of Fig. |I]), then such a solution is allowed only at the 5.3cr CL. In light of the 
strong evidence from SNO that oscillations to a pure sterile state are not responsible for the 
solar anomaly, we hereafter only consider oscillations between active neutrinos. However, 
oscillations to an active-sterile admixture remains viable and are worthy of investigation. 

For a pure active oscillation solution, the NC rate measurement at SNO provides a 
direct determination of the ^B flux produced in the Sun 0. This frees us from relying on 
the Standard Solar Model (SSM) prediction of this flux which has a large uncertainty. That 
is particularly signiflcant because both SK and SNO are mainly sensitive to *B neutrinos. 
These experiments are also sensitive to hep neutrinos, which according to the SSM constitute 
a very tiny fraction compared to the ®B neutrinos. According to Ref. the hep flux is 



in fact the same as that of the SSM with an uncertainty of 20%. To avoid any dependence 
on the SSM ^B flux, previous authors have performed ^B flux free analyses so as to extract 
the oscillation parameters directly from the data pO[. With the recent SNO results, this is 



no longer necessary; the NC rate is itself a B flux measurement and thus a probe of the 
oscillation dynamics. We take this approach. Also, we flx the unoscillated hep flux at the 
SSM value. 

Global analyses of solar data are by now quite standard other than the function 
employed for the statistical treatment |^. We briefly describe the salient features of our 
analysis. We work in the framework of oscillations between two active neutrino flavors and 
focus on the Mikheyev-Smirnov-Wolfenstein (MSW) solutions since the vacuum solution is 



^However, the NC rate derived from SNO's day and night spectra is dependent on the shape of 
the CC energy spectrum, and hence on the oscillation parameters As a result, the CC and NC 
fluxes are strongly statistically correlated 
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tenuous at best; for completeness we also analyze the data in the limit of a pure vacuum 
solution. To evaluate the survival probability of solar neutrinos, we consider neutrino pro- 
duction points that are nonradial and consider the possibihty of double resonances. The 
production point region of the different neutrino fluxes is as given by the SSM. We use 
semi-analytic expressions of the survival probability that have been derived for the almost 
exponential matter density of the Sun ||22|, but numerically integrate the evolution equations 
for the passage of neutrinos through the earth. For the earth-matter density, we make use 



of the Preliminary Reference Earth Model [23|. The time spent by a detector at a particular 



zenith angle is given by the exposure function [^, which determines the extent to which 
earth matter affects the survival probability. 

For the unoscillated neutrino fluxes other than the flux, we adopt the SSM predictions 
from the pp chain and CNO cycle. We treat the flux normalization /3 as a parameter that 
is constrained by the SNO NC measurement. The undistorted spectrum shape of the ^B 
neutrinos is given in Ref. To determine the expected signal at each detector, the fluxes 
are convoluted with the survival probability at the detector, the neutrino cross-sections and 
the detector response functions (for SK and SNO). We use the neutrino-electron elastic 
scattering cross-sections of Ref. ||2^, and the CC and NC cross-sections of neutrinos on 
deuterium of Ref. [E^ . The response functions are given in Ref. [EBl for SK and in Ref. P[ 



for SNO. 

We analyze the event rate from the Homestake experiment [0], the combined rate from 
GALLEX and GNO @,|], the latest SAGE event rate 0, the SK day and night spectra 
corresponding to 1496 days of running 0, and the SNO day and night spectra (which 
include the CC, ES and NC fluxes) J^^- We do not use the SK rate since the SK day and 



night spectra already include the normalization information . 

The statistical significance of an oscillation solution is determined by evaluating a suit- 
ably chosen function. We define a that depends sensitively on the SNO NC flux and 
and is completely independent of the SSM ^B prediction. It is 

x'= E (i?f (/5) - Rm^l)^\RfiP) - R7) 

+ E (^r(/?)-^r)(4i.).?(^f(/?)-^r) 

ij=l,38 

+ E ^r)(4ivo).^-'(^f(/3)-^r)- (9) 

i,j=l,34 

In Eq. d^), i?*^ and R^^^ denote the theoretical and experimental value of the event rate 
or flux measurement (depending on whether i is an experiment or a spectrum bin) nor- 
malized to the expectation for no oscillations. The first term in is the contribution of 
the rate measurements to the analysis. The sum runs from 1 to 3 because the Homestake, 
GALLEX -t-GNO and SAGE rates are included. The 3x3 matrix a\ contains the experimen- 
tal (statistical and systematic) and theoretical uncertainties. This matrix involves strong 
correlations arising from solar model parameters P^]. Note that (j\ does not include a 
theoretical uncertainty for the ®B neutrino flux. 



The second term in Eq. (|^) is the contribution of the distortions of the SK day and 
night spectra and implicitly, that of the SK rate. cr|^ is a 38 x 38 matrix that contains 
the statistical and systematic uncertainties of the 19 + 19 spectral bins. The systematic 
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TABLE II. Best fit values of Am^, tan^ 6 and the ^B normalization /? for each solution from a 
global analysis. The corresponding NC/CC and values are tabulated in the last two columns. 
The number of degrees of freedom is 72. 



Solution 


Am2 (eV2) 




/5 


NC/CC 




LMA 


5.6 X 10-^ 


0.39 


1.09 


3.19 


50.7 


LOW 


1.1 X lO"'^ 


0.46 


1.03 


2.92 


59.9 


SMA 


7.9 X 10-6 


2.0 X 10-3 


1.46 


4.85 


108 


VAC 


1.6 X 10^1° 


0.25 (3.98) 


0.89 


2.46 


76.3 



uncertainties include the energy correlated, energy uncorrelated and energy independent 
contributions. 

The third term encapsulates the contributions of the SNO CC, ES and NC rates and 
distortions of the SNO day and night spectra. The neutron and low energy backgrounds 
are included in i?*''. The flux contribution to Rf^ is multiplied by the normalization 
factor j3 (relative to the central value of the SSM '^B value), which is constrained by the 
oscillation parameter dependent NC and CC flux components of the SNO day and night 
spectra. 

In Fig. 1^, we show the results of an analysis in which Am^, tan^ 9 and (3 have been varied. 
We only show those regions that are allowed at 3cr. Then the Small Mixing Angle (SMA) 
and VAC solutions are absent. We are left with only the LMA and LOW solutions. The 
contours represent the 95.4% C.L. (2(t), 99% and 99.73% C.L. (3ct) allowed regions which 
correspond to Ax^ = 6.17,9.21 and 11.83, respectively. Values of larger than 7r/4 are 
not allowed at the 5cr C.L. The best-fit parameter values from the analysis are presented 
in Table ||. It is significant that the LMA solution is favored over the LOW at the 99% 
C.L. The correlation between the CC and NC rates extractable from the SNO day and night 
spectra goes a long way in choosing the LMA over the LOW. Table ^ shows that the data 
choose large NC/CC ratios in all MSW regions. An NC/CC ratio ~ 3 in the LOW region 
is at the upper end of the range possible in this region. 

The best-fit, minimum and maximum values of the CC day-night asymmetry 0, y4£)Ar(/?), 
at 2cr for the LMA region are displayed in Table fT|. 

We next illustrate how well the best-fits of the two contending solutions, LMA and 
LOW, do in relation to the average survival probabilities of the high energy (*B and hep), 
intermediate energy (''Be, pep, ^^O and ^^N) and low energy {pp) neutrinos extracted from the 



experimental rates. For a description of how these probabilities are obtained see Refs. |13,31 



^The CC day-night asymmetry is defined by 

N-D 

where D and A'' are the total CC fluxes detected during the days and nights, respectively. Discus- 
sions of earth regeneration effects can be found in Ref. pO[|. 
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TABLE III. Best-fit, minimum and maximum values of A]:)n{(3) at 2a for the LMA region. 





Adn{P) 


Am2 (eV2) 




13 


at best-fit 


0.043 


5.6 X 10'^ 


0.39 


1.09 


minimum at 2a 


0. 


1.8 X 10"^ 


0.35 


0.85 


maximum at 2a 


0.122 


2.8 X 10"^ 


0.35 


1.15 



In Fig. ^ we plot each model-independently extracted survival probability at the mean 
energy of the high, intermediate and low energy neutrinos relevant to the experiments. 
The vertical error bars result from the experimental uncertainties in the rate measurements 
and the theoretical uncertainties in the SSM flux predictions. The horizontal error bars 
span the energy ranges of high, intermediate and low energy neutrinos. The solid and 
dashed lines superimposed on the plot are flux-weighted survival probabilities at the SK 
detector corresponding to the best-fit LMA and LOW points of Table [I^, respectively. The 
monoenergetic ^Be and pep fluxes are not included in the averaging. The wiggles in the 
survival probabilities at high energies is a result of earth-matter effects. Aside from the 
averaging, a similar plot was made in Ref. with pre-SNO NC data. 

As a glimpse of the precision in the LMA region that KamLAND data may provide us 
with in three years, we have simulated data (with an antineutrino energy threshold of 3.3 
MeV) at the best-fit LMA point and overlayed the expected 2cr, 99% CL. and 3cr regions 
on the currently allowed LMA region in Fig. ^. 

We conclude that the LMA solution with nonmaximal mixing and a large active neu- 
trino component in the solar flux is favored at the 99% C.L.; the 2a allowed region spans 
2.7 X 10"^ - 1.8 X 10"^ eV^ in Am^ and 0.27 - 0.55 in tan^ 6. KamLAND is an ideal exper- 
iment to precisely measure oscillation parameters in this range. No solution with 6 > tt/A 
is allowed at 5a. Since the mixing is found to be nonmaximal, a determination of the neu- 
trino mass scale via neutrinoless double beta decay is an exciting possibility if neutrinos are 
Majorana A large sterile neutrino component in the solar flux remains viable. 
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2 4 6 

(10^cm"^s"^) 

FIG. 1. Graphical representation of Eq. The diagonal {^nc = ^cc) corresponds to 

sin^ a = or pure sterile oscillations. The dark-shaded band encloses the values of <1> nc a-nd ^cc 
allowed by the SSM at la for sin^ a = 1 or pure active oscillations. The light-shaded band is the 
region allowed by the SSM at la if sin^ q = 0.5. The SNO NC and CC measurements assuming 
an undistorted '^B spectrum (with la uncertainties) are marked with a cross. 
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HOMESTAKE + GALLEX/GNO + SAGE 

+ SK D/N SPECTRA + SNO D/N SPECTRA 

10-3 , , , I , , , , I , , , , I , , , , I , 




0.2 0.4 0.6 0.8 1.0 

FIG. 2. The 2a, 99% C.L. and 3a allowed regions from a fit to the Homestake, GALLEX+GNO 
and SAGE rates, and the SK and SNO day and night spectra. 
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FIG. 3. The flux-weighted survival probabilities of the best-fit LMA (solid) and LOW (dashed) 
solutions from Table |l| in relation to the model-independently extracted values from the data. 
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FIG. 4. Projection of how well KamLAND will determine the oscillation parameters with three 
years of data accumulation assuming an LMA solution. Data were simulated at the best-fit LMA 
parameters. The ellipses are the 2c7, 99% C.L. and 3a KamLAND regions. 
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